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ABSTRACT 

We have analyzed ASCA and ROSAT PSPC spectra and images of the galaxy cluster Abell 2029. The ASCA 
spectra of the cluster indicate that the gas temperature declines with radius. The PSPC image shows that the 
cluster is very regular and smooth. Also, there is no significant evidence for any irregularities in the temperature 
distribution in the cluster, as would be produced by a subcluster merger These results suggest that A2029 is a 
relaxed cluster, and that the gas is in hydrostatic equilibrium. We use the assumption of equilibrium to determine 
the gravitational mass of the cluster as a function of radius. At a radius of 16' (1.92/!5q Mpc; Hg = 50/i5o km 
s~' Mpc"'), the gravitational mass is M,o, = (9.42 ±4.22) x lO'^^/ij^M©, while the mass of gas is Mgas = (2.52 ± 
0.77) X IQ^^k^^^MQ. The gas fraction is found to increase with radius; within a spherical radius of 16', the 
fraction is Mga/M„„ = (0.26 ± 0. U)h-^. 

The iron abundance in the gas is found to be 0.40 ± 0.04 of solar There is no significant evidence for any 
variation in the abundance with position in the cluster. The global X-ray spectra, central X-ray spectra, and 
ROSAT surface brightness all require a cooling flow at the cluster center The global X-ray spectrum implies that 
the total cooling rate is 363'!il1h'^lMQ yr~'. 

The global X-ray spectra are consistent with the Galactic value for the soft X-ray absorption toward the cluster. 
The ROSAT PSPC spectra of the central regions of the cluster are inconsistent with the large value of foreground 
excess absorption found by White et al. (1991) based on the Einstein SSS spectrum. The upper limit on excess 
foreground absorption is 7.3 x 10'^ cm"^. However, the spectra do not rule a significant amount of intrinsic 
absorbing gas located within the cooling flow region. 

Subject headings: cooling flows — galaxies: clusters: general — galaxies: clusters: individual (A2029) — 
galaxies: cD — intergalactic medium — X-rays: galaxies 



1. INTRODUCTION 

Abefl 2029 (A2029) is a richness class II, Bautz-Morgan 
class I cluster of galaxies at a redshift of z = 0.0767 (Abell, 
Corwin, & Olowin 1989). With an X-ray flux of Fx (2-10 keV) 
= 7.52 X 10"" ergs cm"^ s"' (David et al. 1993), it is one of 
the brighter X-ray clusters. This flux corresponds to an X-ray 

"h-l ergs s"' 



luminosity of Lx (2-10 keV) = 2.07 x lO'^^/i^^ ergs s"', and a 



bolometric X-ray luminosity of Lx(bol) = 4.26 x IQ^^h^l ergs 
s"i (Ho = 50/z5() km s"' Mpc"'). 

A2029 is one of the most regular rich clusters known, partic- 
ularly in its X-ray properties (Slezak, Durret, & Gerbal 1994; 
David, Jones, & Forman 1995; Buote & Tsai 1996). As a result, 
it seems likely that the intracluster gas in A2029 is in hydro- 
static equilibrium out to a large radius. Thus, A2029 is an ideal 
candidate for the determination of the total gravitational mass 
of the cluster at large radii using the assumption of hydrostatic 
equilibrium. This requires that the gas temperature and den- 
sity be determined as a function of radius the cluster Previous 
observations using the Einstein Monitor Proportional Counter 
(MPC) found an average gas temperature of kT = 7.8;!;[g keV 



(David et al. 1993). However, this observation does not provide 
any information on the distribution of temperatures. David et al. 
(1995) derived the mass of A2029 from ROSAT Position Sen- 
sitive Proportional Counter (PSPC) observations assuming that 
the gas was isothermal. 

This paper presents new ASCA X-ray spectra of A2029. We 
use these spectra to constrain the temperature variation in the 
cluster. We also analyze archival ROSAT PSPC observations 
of the X-ray image and spectrum of A2029. We use the ASCA 
derived temperature gradients and ROSAT derived X-ray sur- 
face brightness to determine the variation of the total gravi- 
tational mass and gas mass with radius in the cluster As a 
test of the hydrostatic assumption, we search for asymmetries 
in the the temperature distribution which would be indicative 
of recent dynamical activity, such as a subcluster merger We 
also constrain any abundance gradients in the intracluster gas 
in the cluster The use of ASCA spectra to determine the spatial 
variation of the spectral properties of the X-ray emitting gas in 
A2029 requires that the effects of the energy-dependent Point 
Spread Function (PSF) be corrected. This is particularly diffi- 
cult when the X-ray surface brightness is strongly peaked to the 
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cluster center, as is true in A2029. 

A2029 contains a very strong cooling flow, which is centered 
on the large central cD galaxy (UGC 9752 = IC 1101). The 
presence of a cooling flow has been established based both on 
analyses of the X-ray surface brightness profile (Arnaud 1989; 
Sarazin, O'Connell, & McNamara 1992), and on analyses of 
the X-ray spectrum of the central cluster region (White et al. 
1991). The Einstein Imaging Proportional Counter (IPC) obser- 
vations of the surface brightness gave a cooling rate of Mcoi = 
366/!5oM0 yr"^ within a radius of 1.'9 = llGk^Q kpc. (Arnaud 
1989), while the ROSAT High Resolution Imager (HRI) gave 
M„,„, = 370/!J,^Mq yr"^ within 228/z5(i kpc (Sarazin et al. 1992). 
The Einstein Solid State Spectrometer (SSS) spectrum gave a 
cooling rate of Mcooi = 5 1 3;!;247/!5o M0 yr"' within its aperture of 
3' (White et al. 1991). However, A2029 is nearly unique among 
strong, well-studied cooling flows in that it shows no evidence 
for optical emission lines or blue stellar continuum color within 
its central regions (e.g., McNamara & O'Connell 1989). 

White et al. (1991) reanalyzed the Einstein SSS spectra of 
the central regions (3 arcmin radius) of a number of cooling 
flow clusters, and found evidence for large amounts of soft X- 
ray absorption, in excess of that expected from the interstellar 
gas in our Galaxy. In A2029, White et al. found an excess ab- 
sorbing column of ANh = 15 x 10^" cm"^, if the excess absorber 
was modeled as foreground Galactic material. Because the 
Galactic hydrogen column towards A2029 is relatively small, 
Nh = 3.07 X 10^° cm"2 (Stark et al. 1992), this excess absorp- 
tion should have a profound effect on the soft X-ray spectrum 
of the central regions of A2029. 

We extract the central ASCA and ROSAT X-ray spectra of 
A2029 to study the spectral properties of the cooling flow gas. 
In particularly, we use the soft X-ray spectra to constrain the 
amount of absorbing material in excess of that expected from 
our Galaxy. 

In § g, we describe the ASCA and ROSAT observations and 
the treatment of the data. The ROSAT PSPC image of the clus- 
ter field is extracted in § 0; it is needed to analyze the spec- 
tral data and to determine the X-ray surface brightness profile. 
The global X-rayspectrum is analyzed using ASCA and ROSAT 
data in § Q. In § g, we discuss the spatial variation of the X-ray 
spectrum of A2029. Here, we determine the radial tempera- 
ture and abundance gradients, and search for asymmetries in 
the temperature distribution. The central cooling flow spectrum 
is extracted in § ^. Using the X-ray temperature and surface 
brightness profiles, we determine the distribution of the total 
gravitational mass and gas mass in § 0. Finally, our conclusions 
are summarized in § ^. 

Unless otherwise noted, all uncertainties quoted in this paper 
are 90% confidence regions. The abundances given are all rel- 
ative to the "cosmic" values in Anders & Grevesse (1989). All 
distance-dependent values in this paper assume H,, = 50 km s"' 
Mpc"' and qo = 0.0; we have tried to indicate the distance de- 
pendence by parameterizing the Hubble constant as Hg = 5O/150 
km s"' Mpc"^ . 

2. X-RAY OBSERVATIONS 

2.1. ASCA Data 

The cluster was observed with ASCA on 1994 February 19- 
20. The Solid-state Imaging Spectrometer (SIS) detectors op- 
erated with two chips per detector being active (S0C0,1 and 
S1C2,3), and the observation included both bright and faint 
mode SIS data. There were a number of problems with the 



attitude control during a portion of this observation, which re- 
sulted in the original standard attitude file being incorrect. For 
the first 1 180 seconds of the integration, the attitude was about 
20 arcmin off. We discarded this portion of data. The attitude 
settled into an accurate and steady value for the remainder of the 
observation. The Gas Imaging Spectrometer detector 3 (GIS3) 
data was taken during the period when the events were recorded 
with a bit error This error was corrected, although this data 
has fewer spectral channels as a result of this. The data were 
screened using relatively strict versions of the standard clean- 
ing criteria. When these screening operations were complete, 
the exposures were 29,612, 29,710, 34,365 and 34,162 seconds 
for the SISO, SISl, GIS2, and GIS3, respectively. 

The SIS and GIS spectra were corrected for background us- 
ing blank sky observations which were screened in the same 
manner as the data. For the GIS, the background was con- 
structed by weighting blank sky observations with differing 
values of the cut-off rigidity (COR) by the actual distribution 
of COR values during the screened exposure time. For the 
SIS, a long observation with the same minimum cut-off rigidity 
(COR > 8) was used. For the analysis of the spatial variation 
of the X-ray spectrum (§|^), we added a systematic uncertainty 
of 20% in the background to the other uncertainties in deriving 
confidence intervals for parameters. 



2.2. ROSAT PSPC Data 

A2029 was observed with the ROSAT Position Sensitive Pro- 
portional Counter (PSPC) for 3,236 seconds on 1992 January 
24 (observation wp800161; U. Briel, PI.) and for 12,550 sec- 
onds on 1992 August 10-11 (observation rp800249; C. Jones, 
PI.). The first observation was discussed previously in Slezak 
et al. (1994), while the second was presented in David et al. 
(1995). The two observations were combined, and were aligned 
using the positions of two bright point sources. The PSPC 
data was screened for periods of high background based on 
a Master Veto Rate > 170 (Plucinsky et al. 1993), for other 
times of high background based on the light curve of counts 
in the background region used for the extraction of the spec- 
trum (see below), for periods of 15 seconds after switching 
on the high voltage, and for periods with an uncertain aspect 
solution. The analysis of the PSPC was done with a com- 
bination of the IRAF/PROS softwai-e, the XSPEC/XANADU 
software, the XSELECT program, and Snowden's SXRB pack- 
age for the analysis of diffuse X-ray emission (Snowden 1995). 
The screened data contained a total exposure of 12,905 seconds 
and the average Master Veto Rate was 87.05. The data were 
corrected for particle background, scattered solar X-ray back- 
ground, the long term enhancements (LTE) in this background, 
and after-pulses. With these backgrounds removed, the image 
was corrected for exposure and vignetting. 

The ROSAT PSPC spectra and X-ray surface brightness pro- 
files were corrected for background using four regions at radii 
from 24 to 44 arcmin which were free of strong sources or 
obstruction by the instrument ribs. Because of calibration un- 
certainties at the softest and hardest energies, we restricted the 
PSPC spectra to energies of 0.2-2.2 keV (channels 6-32). 



2.3. Spectral Modeling 

We fit the observed spectra with thermal emission models 
for a low density plasma in coUisional ionization equilibrium. 
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Where possible, we used models based on both the Raymond- 
Smith code (Raymond & Smith 1977) and the MEKAL code 
(e.g., Kaastra et al. 1996). The MEKAL code includes a better 
treatment of the iron L lines near 1 keV, based on measurement 
and calculations of the atomic parameters by Liedahl and his 
colleagues (e.g., Liedahl et al. 1995). We did not fit MEKAL 
models during the spatial analysis of the spectra in § ||, as the 
MEKAL spectra were not available in a form which could be 
used with the programs employed to correct for the energy- 
dependent Point-Spread-Function (PSF). However, at higher 
temperatures, we found that the Raymond-Smith and MEKAL 
spectra gave almost exactly the same results. Larger differences 
are found between the Raymond-Smith and MEKAL spectra at 
lower temperatures, and these become more important when 
the spectrum includes the cooling flow regions at the center of 
the spectra. However, we still found that the spectral parame- 
ters from fitting Raymond-Smith and MEKAL models were in 
good agreement. Because the MEKAL model generally gave 
slightly better fits, we report the results using this model in §§ Q 
and|. 

We considered two simple models for the thermal distribu- 
tion of the gas in each region where the spectrum was extracted. 
The first model assumed thermal emission at a single temper- 
ature; we refer to this as a "single temperature" model. In Ta- 
ble |l| these models labeled as "IT" This model is characterized 
by two spectral parameters, the temperature of the gas and the 
ratio of the abundances in the gas to the solar values, and by 
the overall emission integral of the gas which normalizes the 
spectra. 

As discussed above (§ |l]), there is considerable previous ev- 
idence from X-ray spectra and images for a cooling flow at 
the center of A2029. Thus, we also considered a simple spec- 
tral model for ambient gas and a cooling flow. The gas in the 
cooling flow was assumed to cool from the temperature of the 
ambient gas to very low temperatures due to its own radiation. 
The abundances in the cooling flow gas were taken to be those 
in the ambient gas, as expected if the cooling gas originated 
as intracluster gas. The cooling flow gas was assumed to cool 
at constant pressure; more detailed spectral models have indi- 
cated that this is a reasonably accurate approximation (Wise & 
Sarazin 1993). The X-ray spectrum of this cooling flow com- 
ponent is given by 
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where E = hv is, the photon energy, LecIE is the luminosity 
emitted in the photon energy range E -^ E + dE, /i « 0.61 is 
the mean mass per particle in terms of the proton mass nip, T is 
the ambient gas temperature in the cluster from which the gas is 
assumed to cool, T' is its current temperature, the emissivity of 
the gas per unit volume and photon energy is p^Ae{T), p is the 
mass density of the gas, and the total emissivity (cooling rate) 
of the gas is p^A,o,(T) (e.g.. Wise & Sarazin 1993). We wifl 
refer to this model as the "single temperature plus cooling flow 
model." In Table nl these models labeled as "CF." This model is 
characterized by a single additional parameter, the cooling rate 

Mcool- 

The global cluster X-ray spectrum (m and the spectrum of the 
inner cooling flow region (§0) were fit to these models using 
XSPEC version 9. Because this version of XSPEC did not con- 
tain a model for the cooling flow spectrum using the MEKAL 
emission code, we constructed such a model. The analysis of 



the spatial variation of the spectrum, correcting for the effects 
of the energy-dependent PSF of ASC A, was done with the pro- 
gram of Markevitch et al. (1996). 

3. ROSAr PSPC IMAGE 

The ROSAT PSPC data in flie hard band of 0.5 to 2.0 keV 
(PI channels 52 to 201) were used to construct an image of 
the region of the ASCA observation. This hard band is useful 
because it reduces the background in the image without much 
loss of counts from the hot intracluster gas. This image was 
corrected for particle and solar X-ray background, exposure, 
and vignetting as discussed above. The resulting image was 
adaptively smoothed with a variable smoothing kernel so that 
each smoothing beam had a signal to noise ratio of at least five 
(Huang & Sarazin 1996). Contours of the X-ray surface bright- 
ness for the central portion of the resulting image are shown in 
Figure |l|. 
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Fig. L — Contours of the ROSAT PSPC image in tlie hard band (0.5 to 2.0 
keV) of the A2029 cluster field, corrected for particle and solar X-ray back- 
ground and for exposure. The image was adaptively smoothed to a signal to 
noise ratio of five per smoothing beam. The contours are logaiithmically spaced 
with four contours per dex, and the lowest contour coiresponds to 0.001 counts 
s"' arcmin"-. The coordinates are J2000. Two unrelated point sources located 
approximately 4 and 9 arcmin from the cluster center are labeled as "A" and 
"B," respectively. The two dashed squares give the approximate fields of view 
of the two chips of the SISO detector on ASCA. 



The ROSAT PSPC image of A2029 contains a large number 
of presumably unrelated point sources. There are two moder- 
ately strong point sources located approximately 4 arcmin (Src. 
A) and 9 arcmin (Src. B) from the center of the cluster Src. 
A is known to be a background active galaxy nucleus (AGN; 
Bregman, private communication). Src. B has not been identi- 
fied to our knowledge. These sources are labeled on Figure nl 
and are discussed further in § pi where they are included in the 
determination of the variation in the X-ray spectrum of A2029. 

When these sources are removed, the contours of the X- 
ray surface brightness of A2029 in Figure H appear ellipti- 
cal and quite regular This agrees with the previous conclu- 
sions reached by Slezak et al. (1994) using wavelet techniques 
applied to the shorter ROSAT PSPC exposure, and by Buote 
& Tsai (1996) using the power ratio technique applied to the 
longer ROSAT PSPC exposure. There is general agreement that 
A2029 is regular and that it lacks any strong substructure. This 
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is also consistent with the tendency for clusters with large cool- 
ing flows to be regular (e.g., Buote & Tsai 1996). The regularity 
of A2029 makes it a particularly attractive candidate for apply- 
ing the hydrostatic method to determine the mass distribution. 
At the faintest levels, the ROSAT PSPC observation shows 
X-ray emission from A2029 extending out to a radius of ap- 
proximately 20 arcmin from the cluster center. In § |[ we will 
analyze the X-ray spectrum of A2029 out to a radius of 16 ar- 
cmin. All of this region is included in the fields of view of both 
of the GIS detectors. However, the SIS detectors only cover a 
fraction of the cluster emission. In Figure [|, the dashed squares 
give the approximate fields of view of the two chips used on the 
SISO detector (SOCO and SOCl). The corresponding regions for 
the SISl detector chips (S1C2 and S1C3) are similar and have 
about the same orientation, but are displaced about 1 arcmin to 
the west. 



4. GLOBAL X-RAY SPECTRUM 

The ASCA and ROSAT observations were used to determine 
the integrated spectrum of the cluster For this purpose, we de- 
termined the spectrum within 16' (1. 92/150 Mpc) of the center 
of the cluster This region is shown as the outermost dashed cir- 
cle in Figure |[ This region extends nearly to the bright ring of 
enhanced background near the edge of each of the ASCA GIS 
detectors, and lies well within the ribs on ROSAT PSPC detec- 
tor Since the ASCA SIS spectra cover only a portion of this 
region (Fig. [l|), we do not include the SIS data in this deter- 
mination of the global spectrum. Later, we will find that the 
cluster has significant spatial variations in its X-ray spectrum, 
due to a cooling flow at the c lust er center (§ ^ and due to a 
radial temperature gradient (§ 5.1). Thus, it would be difficult 
to directly compare the SIS spectra with those from the GIS 
and PSPC, which cover the same larger region of the cluster 
The spectra were binned to insure that there were at least 20 
counts per spectral channel, so that the x^ test gave reliable re- 
sults. We excluded appropriately sized regions around each of 
the two point sources labeled A and B in Figure y 

The results of fitting "single temperature" and "single tem- 
perature plus cooling flow" models to the observed spectral are 
summarized in Table |l| The first column enumerates the rows 
for reference in the text. The second column indicates the re- 
gion of A2029 being considered, and the third column lists the 
instrument(s) used. The fourth column gives the model being fit 
to the data. Columns 5-8 give the galactic hydrogen absorbing 
column A^H, the ambient gas temperature kT, the heavy element 
abundance, and the cooling rate (the latter is only given for the 
CF models). Column 9 gives the resulting x^ and number of 
degrees of freedom (d.o.f.) for the best-fit model. The final 
column gives the probability that the value of x^ would be this 
large if the data were drawn from the best-fit model with only 
statistical uncertainties. Because it is likely that there are sys- 
tematic errors in the responses of the instruments, these num- 
bers may not be accurate. The uncertainties in the parameters 
give the 90% confidence region. 

Because of their nearly disjoint spectral bands, we find that 
the ROSAT PSPC (0.2-2 keV) and ASCA GIS's (0.5-11 keV) 
provide nearly orthogonal but complementary information of 
the spectrum of A2029. We first discuss the results from ROSAT 
and ASCA separately, and then jointly fit all of the spectra. 



4. 1 . ROSAT PSPC Global X-ray Spectrum 



The ROSAT PSPC spectrum of the inner 16' of A2029 was 
extracted and corrected for background as discussed above and 
in § 2.2. First, a single temperature model was fit to this spec- 
trum; the results are shown in the first row of Table [l|. The fit 
was not acceptable, with a x^ = 34.05 for 23 d.o.f. This unac- 
ceptable fit gave a well-defined value for the absorbing column, 
but fairly uncertain values for the temperature and abundance 
of the gas. For the high ambient gas temperature in A2029, 
the low energy band of ROSAT makes it relatively insensitive 
to changes in the temperature. Moreover, at these high temper- 
atures, the strongest spectral feature due to heavy elements is 
the 7 keV iron K line complex, which is beyond the band of 
ROSAT. On the other hand, the low energy band of the PSPC 
does make it very sensitive to the Galactic absorbing column. 
The absorbing column required by the PSPC spectrum is con- 
sistent with the Galactic value of A^h = 3.07 x 10^° cm"^ (Stark 
et al. 1992). However, the temperature is lower than that given 
previously by the Einstein MPC (kT = 7.8+| Jj keV; David et 
al. 1993) or that determined below from ASCA (kT = 8.69!j]ji] 
keV; § 4.2 ). An examination of the residuals for this single tem- 
perature fit suggested that the poor fit was due to additional soft 
X-ray emission. 

The poor fit of a single temperature spectrum to the PSPC 
data and the previous evidence for a cooling flow at the center 
of the cluster suggest that a cooling flow component be added to 
the model spectrum. The best-fit single temperature plus cool- 
ing flow model is shown in the second row of Table |l[ This 
provided a considerably improved fit (a reduction of Ax^ w 10 
for one additional d.o.f.). In this model, 53% of the counts in 
the PSPC spectrum come from the cooling flow component. 
Because of the soft band of the PSPC, the temperature of the 
gas is very poorly defined in this model. At the 90% confidence 
limit, only a lower limit on the temperature is determined. As 
noted before, varying the temperature has little effect on the 
soft X-ray spectrum in the single temperature model, and even 
less in the cooling flow model (eq. |l]). The reason for the latter 
is that most of the softer emission in the cooling flow model 
comes from gas at which has cooled to lower temperature, so 
that the integral in equation (|l]) is nearly independent of the up- 
per limit temperature T for E -^kT. 

The abundance of heavy elements is also rather poorly deter- 
mined by the single temperature plus cooling flow fit to spec- 
trum, and only gives a lower limit. As noted above, the Fe K 
lines cannot be observed with ROSAT to constrain the abun- 
dances. The cooling flow component does produce strong line 
features within the PSPC band due to heavy elements, includ- 
ing the L lines of iron and the K lines of C, N, O, Ne, and Mg. 
However, the strengths of these lines are not strongly dependent 
on the abundance of the heavy elements, at least if they are all 
assumed to vary together and the abundance is not too low. The 
reason is that, at temperatures kT ^ I keV where these lines 
are strongest, heavy element line emission dominates the cool- 
ing of the gas. Thus, if the abundances are increased, the line 
emission rate is increased but the cooling time is decreased by 
nearly the same factor In equation (p, increasing the abun- 
dance increases both A£(r') in the numerator and A,o, in the 
denominator by approximately the same amount. 

On the other hand, the ROSAT PSPC spectrum does give an 
excellent determination of the absorbing column, and a good 
determination of the cooling rate Mc„oi- As noted above, the ab- 
sorbing column is reasonably close to the Galactic value from 
Stark et al. (1992). The cooUng rate also agrees well with previ- 
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Table 1 



ASCA AND ROSAT X-RAY SPECTRAL Fits 



Row 


Region 


Instr." 


Model" 


Nh 


kT 


Abund. 


Ma,„ihio 


x'/d.o.f. 


Prob. 










(10-" cm^^) 


(keV) 


(0) 


{Me yr-') 




(%) 


1 


Whole Cluster 


PSPC 


IT 


T 9, +0.22 
J-^i-071 


D-13_o.82 


0.99^:^1 


— 


34.0y23 = 1.48 


6.44 


2 


(r<16') 


PSPC 


CF 


9 S 1+0-20 
^•01-0.16 


>6.23 


>0.34 


377!?^ 


24.47/22=1.11 


32.31 


3 




GIS 


IT 


<4.43 


8.69tiJl 


n -37+0-04 

"•J '-0.03 


— 


788.0^727=1.08 


5.77 


4 




GIS 


CF 


<6.04 


9.54!!:« 


0.40 ±0.04 


<609 


786.0^726=1.08 


6.05 


5 




PSPCGIS 


IT 


3.21±0.11 


8.06±0.21 


0.38 ±0.03 


— 


860.1/753 = 1.14 


0.64 


6 




PSPCGIS 


CF 


3.29±0.12 


9.35tl? 


0.40 ±0.04 


363_9s 


821.9^752=1.09 3.87 




7 


0' <r< 1.5' 


ASCA 


CF 


3.07 


9A2-\f, 


0.40 


216_i2o 


172.61/197 = 0.88 




8 


1.5' <r<5' 


ASCA 


IT 


3.07 


SM^Jit 


0.40 




— 




9 


5' <r< 10' 


ASCA 


IT 


3.07 


6.35lil 


0.40 




— 




10 


10'<r< 16' 


ASCA 


IT 


3.07 


6.295:§? 


0.40 




— 




11 


-20° <PA< 70° 


ASCA 


IT 


3.07 


8.85l;Ji 


0.40 




173.70^212 = 0.82 




12 


70° <PA< 160° 


ASCA 


IT 


3.07 


5.65!^:?^ 


0.40 




— 




13 


160° < PA < 250° 


ASCA 


IT 


3.07 


n TO+4.10 
'•^0-2.70 


0.40 




— 




14 


250° <PA< 340° 


ASCA 


IT 


3.07 


n.orZl" 


0.40 




— 




15 


Center 


GIS 


IT 


T -,£- + 1.41 


i.eo^i] 


0.38 ±0.05 





624.1/557 = 1.12 


2.53 


16 


(r < 3') 


SIS 


IT 


6.64 ±0.38 


f. c^+0.18 

6.55_o.i7 


0.41 ±0.04 


— 


1018.0^816=1.25 


<0.01 


17 




PSPC 


IT 


3.25^:1^ 


4.75!°- 


0.37!°;M 
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27.1/23 = 1.181 


25.17 
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GIS 


CF 


11.45!!;°? 


9.611:?1 


0.38 ±0.05 
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618.2/556=1.11 


3.44 
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SIS 


CF 


V.62!«i° 


7 61+0-45 
'•Oi-0.38 


0.41 ±0.04 


200!ii 


975.30^815 = 1.20 


<0.01 
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PSPC 


CF 


2.89 ±0.23 


>6.36 


>0.30 


329^^6 


20.1/22 = 0.91 


57.68 


21 




All 


CF 


3.29 ±0.09 


1 co+0.28 
'•JO-0.16 


0.43 ±0.03 


<32 


1959.yi401 = 1.40 


<0.01 



"The instruments are: PSPC: R05ATPSPC; GIS: ASCA GIS2 and GIS3: SIS: ASCA SISO and SISl: ASCA: ASCA GIS2, GIS3, SISO, & SISl; Ah: PSPC, SIS, & GIS. 
*1T is the single temperature model, while CF is the single temperature plus cooling flow model. 



ous determinations; we discuss this in more detail below (§ 4.3 ) 



4.2. ASCA GIS Global X-ray Spectrum 

The ASCA GIS2 and GIS3 spectra in the photon energy range 
of 0.5-1 1 keV were determined for the inner 16' of the cluster 

We didn't use the ASCA SISO or SISl spectra, because these 
instruments did not cover the entire region of the cluster (Fig- 
ure nl). Later, we will find that the cluster has significant spatial 
variations in its X-ray spectrum, due to a cooling flow at the 
cluster center (§ ^ and due to a radial temperature gradient 
(§ 5.1). Thus, it would be difficult to directly compare the SIS 
spectra with those from the GIS and PSPC, which cover the 
same larger region of the cluster 

The results of fits to these spectra are summarized in rows 
3 and 4 of Table n^. Because of the hard energy band of the 
GIS spectra, these observations have strengths and weaknesses 
which are the opposites of those of the ROSAT PSPC. The GIS 
provides accurate measures of the ambient gas temperature and 
the iron abundance, but do not usefully constrain the soft X-ray 
absorption or the cooling rate. 

These two spectra were fit first with the single temperature 
model. This gave an accurate determination of the temperature 
and of the iron abundance, but only an upper limit on the ab- 
sorbing column. The fit was adequate, with a x^ per d.o.f. of 
1.08. Because of the evidence for a cooHng flow in A2029, we 



also tried a fit which included a cooling flow. This provided 
only a marginal improvement in the fit (Ax^ = 2 for one extra 
fitting parameter), and accordingly, only an upper limit on the 
cooling rate. 

4.3. Combined ASCA-ROSAT Global Spectral Fits 

Because of the complementary strengths and weaknesses of 
the the ASCA and ROSAT spectra, and the fact that the best- 
fit individual spectra gave consistent results (Table M, we de- 
termined the best-fit spectrum for the combined data sets of 
ROSAT PSPC and ASCA GIS2 and GIS3 spectra. We allowed 
for a variation in the overall calibration between ASCA and 
ROSAT, but required that the shape of the spectrum agree. The 
results of fitting single temperature and single temperature plus 
cooling flow models to the spectra are shown in rows 5 and 6 of 
Table ll|. The abundances in these fits are essentially the same 
as those given by the GIS spectra alone, while the temperatures 
are slightly lower The absorbing columns are similar to those 
from the PSPC alone, but with smaller uncertainties. The cool- 
ing rate in the CF model is basically that same as that given by 
the PSPC alone. 

The addition of a cooling flow improves the fit significantly 
(Ax^ = 38 for one extra fitting parameter). This is due almost 
entirely to the effect of this component on the fit to the PSPC 
spectrum. Because the single temperature plus cooling flow 
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model provides fairly good fit to both the ROSAT and ASCA 
data, we will adopt the spectrum in this fit as our fiducial global 
cluster spectrum for A2029. The upper panel in Figure g shows 
the observed global ASCA GIS2 and GIS3 spectra and this best- 
fit model; the lower panel shows the residuals to the fit. Figure || 
shows the ROSAT PSPC spectrum, model, and residuals. 

The hydrogen absorbing column required by this fit is in rea- 
sonable agreement with the Galactic value of A^h = 3.07 x 10^" 
cm"^ from Stark et al. (1992). The temperature in the single 
temperature plus cooling flow model is higher than the value of 
kT = 7.8;|;^ keV found by David et al. (1993) with the Einstein 
MPC, although the error bars overlap. However, David et al. 
fit a single temperature model to their spectrum. If one com- 
pares to the single ternperature joint fit to the ASCA and ROSAT 
data (row 5 in Table [I]), the agreement is quite good. The to- 
tal cooling rate is in remarkably good (and probably fortuitous) 
agreement with the values of M, - t'^'^i^-^ n/f ^----i 



^cool 



366/!5oM0 yr ' derived 



with the Einstein IPC (Arnaud 1989), and M„,„/ = 37O/i;oM0 
yr"' derived from the ROSAT HRl (Sarazin et al. 1992). Both 
of these determinations were based on the surface brightness 
profile of the cluster, rather than the X-ray spectrum. The cool- 
ing rate was also derived using the X-ray spectrum from the 
Einstein SSS by White et al. (1991). Their value of M^ooi = 
^^^-liVho^Q y' is consistent with ours within their very 
large uncertainties. 

The X-ray flux in these models is Fx (2-10 keV) = (1.00 ± 
0.01) X 10""^ ergs cm"^ s"'. The uncertainty shown is statisti- 
cal, and is dominated by the good statistics in the ASCA GIS 
spectra. However, the systematic uncertainty in the absolute 
calibration of these ASCA spectra is much larger This flux is 
considerably larger than that found by David et al. (1993), Fx 
(2-10 keV) = 0.75 x 10"'" ergs cm"^ s"', using the Einstein 
MPC, which had a larger field of view than ASCA. 

5. SPATIAL VARIATION OF THE X-RAY SPECTRUM 

The ASCA observations were used to study the spatial vari- 
ation in the X-ray spectrum, and derive the radial variation in 
the gas temperature and in the abundances. Both SIS and GIS 
data were used. The determination of the spatial variation in 
the spectrum is complicated by the broad, energy-dependent 
Point-Spread-Function (PSF) of the ASCA mirrors. We used 
the technique of Markevitch et al. (1996) to correct for the re- 
distribution of X-ray photons produced by the ASCA PSF. A set 
of "model" regions on the sky were selected, and we assumed 
that the shape of the X-ray spectrum was constant within each 
of these regions. A corresponding set of "image" regions were 
defined in the focal plane of ASCA. Generally, the image re- 
gions and model regions were identical, except that the outer- 
most model region was larger than the outermost image region. 
This was done to avoid the unphysical effects of assuming that 
the cluster emission ended abruptly at the outer edge of the last 
image region. Also, the image regions were restricted to the 
area within 18' of the center of each of the GIS detectors. 

The observed ASCA X-ray spectra were accumulated within 
each image region. For a given set of spectral parameters for the 
model regions, we determined the predicted spectra in the im- 
age regions, including the detector response and the redistribu- 
tion of photons among image regions by the energy-dependent 
ASCA PSF. The parameter of the spectra in the model regions 
were varied until the best-fit model (minimum x^) for the ob- 
served spectra in the image regions was found. An annealing 
technique (Press et al. 1992) was used to avoid false minima. 



We were careful to use spectral bins which were large enough 
that the Gaussian approximation to Poisson uncertainties was 
valid, so that the x^ test could be applied. 

The energy-dependent ASCA PSF was approximated by in- 
terpolation between GIS images of Cyg X-1 taken at various lo- 
cations in the focal plane (Takahashi 1995). The PSF is poorly 
known at photon energies below 2 keV. For photon energies be- 
tween 2.0 and 2.5 keV, the effective area of ASCA is affected 
by a poorly calibrated absorption edge from gold in the mir- 
rors. Thus, we have chosen to fit only photon energies above 
2.5 keV in the spatial analysis. The lower energy cutoff should 
not affect the determinations of the high temperatures found for 
the gas. Excluding lower energies does reduce the sensitivity 
of ASCA to the central cooling flow (§ ^, but also reduces the 
calibration problems encountered with ASCA at low energies 
(mainly below 1 keV; § 5.1). A relative systematic PSF un- 



certainty (5% for circularly symmetric image regions, 15% for 
other region shapes) was added to the statistical uncertainties. 

The instrumental contributions to the PSF are different for 
the SIS and GIS. For modeling of the SIS data, the GIS Cyg X- 
1 calibration images were corrected for the energy dependence 
of the intrinsic GIS detector blurring by additional smoothing 
resulting in a final constant resolution (Gaussian a = 0.'5). The 
cluster SIS data also were smoothed to the same resolution. 

Because of its superior spatial resolution, the ROSAT PSPC 
image (Fig. |l]) was used to determine the distribution of pro- 
jected emission measure within the cluster. Errors due to pho- 
ton statistics and the background correction in the ROSATPSPC 
image where included in the spectral analysis. At each step 
in the iteration to find the best-fit ASCA spectra, the spec- 
trum in each model region was used to determine the emis- 
sion within the band of the ROSAT PSPC image including the 
ROSAT PSPC response, and the surface brightness in the PSPC 
image was converted into a distribution of emission measure on 
the sky. Because there is some uncertainty in the overall ROSAT 
to ASCA calibration, we allow this to vary, and only use the 
ROSAT PSPC image to fix the relative emission measure varia- 
tion across the sky. This emission measure distribution and the 
spectral parameters were used to determine the predicted A5CA 
spectra within the image regions, including the redistribution 
by the PSF. Because a small miss-alignment of the ROSAT and 
ASCA images can affect the spectral fit, it is important to align 
the ROSAT and ASCA data accurately. We did this, but we also 
vary the offset between the two satellite data sets, and include 
an estimate of the resulting uncertainties in the alignment on 
our spectral fits. 

The use of the ROSAT PSPC image to constrain the emission 
measure distribution in the ASCA only works well if the same 
spectral component provides the bulk of the emission in these 
two disjoint spectral bands (0.5-2.0 keV for the ROSAT image, 
and 2.5-1 1 keV for the ASCA data). This is probably a rea- 
sonable assumption in the outer parts of A2029, but not in the 
inner cooling flow region. There, most of the ROSAT photons 
come from cooling gas, while much of the ASCA spectrum is 
contributed by hotter, ambient gas. To avoid errors due to this, 
we free the normalization of the central ASCA spectrum relative 
to the outer regions during the spectral fitting. 

Two bright X-ray point sources are seen on the ROSATPSPC 
image (Fig. [l|) at about 4 and 9 arcmin from the cluster center 
A region around each of these sources was excluded from the 
ASCA cluster spectra. However, because of the broad ASCA 
PSF, some emission from each of these point sources may be 
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Fig. 2. — The global ASCA GIS2 and GIS3 X-ray spectra for A2029 (r < 16') as a function of the measured photon energy are 
shown is the upper panel. The crosses give the data points with l-cr error bars, while the histogram is the best single temperature 
plus cooling flow model fit to the PSPC, GIS2, and GIS3 spectra (row 6 in Table nh. The width of the data points or histogram steps 
is the width of the energy channels used to accumulate the data. The upper spectrum is for the GIS3, while the lower spectrum is fo r 
GIS2. The GIS3 has larger energy bins because the bit error in memory during the observation reduced the energy resolution (§ 2.1). 
The lower panel gives the residuals to the fit (in counts/sec/keV). 



included in the adjoining cluster image regions. Thus, we ac- 
cumulate the ASCA X-ray spectrum from the image regions 
around each of these sources, and simultaneously fit a power- 
law spectrum to each of these sources. (Thermal spectra were 
found to provide poorer fits than power-law spectra.) We found 
that contamination from these sources only has a small effect 
on the derived spectral properties of the cluster regions except 
for the outermost cluster regions. Because these sources may 
have varied between the time of the ROSAT and ASCA obser- 
vations, we allow the overall normalization of the flux of these 
sources to vary. In fact, the ASCA and ROSAT fluxes of source 
A (Fig. IT]) agree within the uncertainties when corrected for the 
different spectral bands, but source B is about 3 times fainter 
in the ASCA observation than in the ROSAT data. The pho- 
ton power-law indices of sources A and B were -1.49;!;n^T and 



-1.1 l;!;o34, respectively, in the best-fit models. 



-0.42 



The uncertainties in derived spectral quantities are 90% con- 
fidence regions from 200 Monte Carlo simulations of fitting the 
data with random statistical and systematic uncertainties added. 



5.1. Radial Temperature Gradient 

The radial temperature variation in A2029 was determined 
by accumulating the ASCA SIS and GIS spectra in four adjoin- 
ing annuli with boundary radii of 0, 1.5, 5, 10, and 16 arcmin. 
The outer radius of the outermost model annulus was taken to 
be 20 arcmin. The spectra were grouped into eight spectral 
channels with boundary photon energies of 2.5, 3.0, 3.5, 4.0, 



5.0, 6.0, 7.0, 8.5, 11.0 keV in the GIS detectors. In the SIS 
detectors, the last two of these channels were combined into 
a single channel from 7 to 1 1 keV. For the purposes of deter- 
mining the temperature variation, the abundances of the heavy 
elements were fixed at 0.40 of the solar value, which was the 
best-fit abundance for the entire cluster spectrum (§ 0). 

We initially tried to fit the spectra assuming a single temper- 
ature model for the spectrum in all of the cluster annuli. How- 
ever, this did not provide a very good fit to the central circular 
region. Even when we freed the normalization of this region 
(i.e., no longer set the normalization of this region relative to 
the others from the ROSAT PSPC image), the fit was only im- 
proved very slightly. Examination of the residuals indicated 
that the model was unable to provide enough soft flux in the 
central region. Previous X-ray spectral (e.g.. White et al. 1991) 
and spatial observations (e.g., Sarazin et al. 1992), and the cen- 
tral ROSAT PSPC and ASCA spectra presented in this paper 
(§ 0), all indicate that A2029 contains a cooling flow in its cen- 
tral regions. Therefore, we added a cooling flow component to 
the central circle region in all of our spatially -resolved spectral 
fits. The gas was assumed to cool from the ambient temperature 
in this region (the temperature of the Raymond-Smith compo- 
nent), and have the same abundances as the intracluster gas. 
The addition of this cooling flow component lowered \^ by 32 
for only one additional parameter. This model did provide an 
acceptable fit. 

The best-fit spectral parameters for this model for the cluster 
annuli are shown in rows 7-10 of Table [1] Note that the best-fit 
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Fig. 3.— The global ROSAT PSPC X-ray spectrum for A2029 (r < 16'^ 
the upper panel. The crosses give the data points with l-a error bars, whil^the 
flow model fit to the PSPC, GIS2, and GIS3 specti-a. (row 6 in Table [I]). T^ 
of the energy channels used to accumulate the data. The lower panel gives 
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as a function of the measured photon energy is showt 

histogram is the best single temperature plus cooling 

width of the data points or histogram steps is the wi<lth 

the residuals to the fit (in counts/sec/keV). 
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value of x^/d.o.f. is significantly less than unity. This is due 
to the inclusion of systematic uncertainties in the evaluation of 
X^, as discussed above. As a result of this, we do not give the 
probabilities associated with these fits. 

The best-fit temperatures in annuli are also plotted in Fig- 
ure 0. There is some evidence for a decrease in the temperatures 
with radius. This trend was found in all of the fits to the ASCA 
spectra, including the fits which did not correct for the two 
unrelated sources or include the effects of the central cooling 
flow. To assess the significance of the temperature gradient, the 
spectra in annuli were fit assuming that the temperature was the 
same at all radii. The horizontal dashed line in Figure gives 
this best-fit isothermal model for all of the spatially-resolved 
spectra {kT = 8.6 keV). Note that this temperature is slightly 
different than the best-fit temperature for the whole cluster (also 
including the effects of a cooling flow), ^r = 9.35 keV, for sev- 
eral reasons. First, This was a much worse fit to the spectra. 
Based on the F-test and the simulations of the uncertainties in 
the fits, this isothermal fit can be rejected at the > 96% confi- 
dence level. To parameterize the temperature gradient, we also 
fit the spectra in annuli assuming that the temperature was a lin- 
ear function of the radius. The slanting dashed line in Figure 
is the best-fit linear function of radius. 
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Fig. 4 

of radius as derived from the ASCA spectra. The vertical error bars give the 
90% confidence region for the temperature. The horizontal error bars show the 
radii of the annuli over which the spectra were accumulated. The temperature is 
plotted at the median radius of the X-ray emission in that annulus. The horizon- 
tal dashed line gives the best-fit single temperature fit to the individual annular 
spectra. The sloping dashed line is the best-fit linear function of the radius 
(eq. Q). The spectral model used for the central bin includes a cooHng flow, and 
the temperature plotted is the upper limit temperature, which we assume is that 
of the ambient gas in this region. 



Fits to the GIS data alone showed the same declining trend 
of temperatures with radius and consistent values for the tem- 
peratures. In fits with the SIS data alone, the temperatures in 
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It is an improvement over the isothermal temperature distribu- 
tion, with a reduction in y^ of 9 for only one additional param- 



tion, with a reduction in x 
eter. 
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the inner two annuli were similar to those in the combined fit. 
The temperatures in the outer two annuli were very poorly de- 
termined, due to the small portion of these regions included on 
the SIS chips (Figure |l]). However, the temperatures were con- 
sistent within the uncertainties, which were very large in the 
outer two annuli. 



5.2. Radial Abundance Variation 

We also derived the abundances in the same radial rings 
which were used in the previous section. The values of the 
derived abundances are given in Figure ^ Unfortunately, the 
uncertainties in the individual abundances are quite large, and 
only upper limits are found for the outer two annuli. The is 
not really due to the statistical uncertainties on individual spec- 
tra. Instead, it is largely a consequence of the energy dependent 
PSF of ASCA. The PSF is particularly wide at high energies, 
so the Fe K line complex, which is the primary spectral feature 
determining the abundance, can be strongly affected. When 
combined with the very strong central surface brightness of the 
Fe K line complex in this cooling flow cluster, this implies that 
there is a very large correction for scattering of the central emis- 
sion in the Fe K line spectra of the outer regions. This results 
in very large uncertainties for the outer abundances, and causes 
the uncertainties in the outer abundances to be anticorrelated 
with the uncertainties in the central abundance. 




arcmin) 



Fig. 5. — The iron abundance of the intracluster gas in A2029 as a function 
of radius as deiived from the ASCA spectra. The vertical enor bars give the 
90% confidence region for the abundances. The horizontal error bars show the 
radii of the annulus over which the spectra were accumulated. The abundance is 
plotted at the median radius of the X-ray emission in that annulus. The dashed 
horizontal line is the best-fit single abundance value for the cluster There is no 
statistically significant variation of the abundance with radius. 

Figure Bl shows some slight evidence for either a decline 
in the iron abundance with radius or a higher abundance in 
the center cooling flow region than in the rest of the cluster. 
However, the difference between the spectral fit with individual 
abundances and a fixed overall abundance is not statistically 
significant (a decrease in x^ of 3.3 for three additional param- 
eters). One finds that the hypothesis that the abundance is con- 
stant can only be rejected at the 51% confidence level (i.e., not 
at all). For this reason, we have assumed that the abundance is 
constant at the value set by the global cluster spectrum in the 
determinations of the cluster temperature variation. 



5.3. Azimuthal Temperature Variation 

We also determined the azimuthal variation in the X-ray tem- 
perature in A2029. As shown in Figure U the X-ray emission 
in A2029 is elongated at a position angle (PA) of approximately 
25° (measured from the north to the east). We divided the re- 
gion of the cluster between radii of 3' and 16' into four annular 
sections, each 90° wide, with the first centered on PA = 25°. 
The temperature was determined independently for the central 
3 arcmin radius circle and for each of these annular sections. 
These regions are shown in Figure ^ As before, regions were 
used to isolate the two point sources A and B. The abundance 
was fixed at 0.40 of solar, the best-fit value for the cluster spec- 
trum derived above. The relative systematic PSF uncertainty 
was increased to 15% because of the shape of the regions used. 
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Fig. 6. — The annular sections used to derive the azimuthal variation of the 
temperature in A2029. The contours give the ROSAT PSPC surface brightness 
as in Figure M. The outer dashed circle has a radius of 16 arcmin. Only regions 
interior to this circle were used to derive the global or spatially resolved X-ray 
spectra. The inner dashed circle has a radius of 3 arcmin; the spectrum interior 
to this regions was used to determine the properties of the cooling flow in the 
cluster (§ rl). The radial dashed lines delimit the 4 regions used to determine 
the azimuthal temperature variation (§pj). 



The resulting temperatures for the four regions are listed in 
rows 11-14 of Table [u, and are shown as a function of the po- 
sition angle in Figure^ There is some evidence that the gas is 
hotter to the northwest, and cooler to the southeast. We also 
fit a single temperature to all four sectors; this is shown as 
a dashed horizontal line in Figure 0. Formally, the measured 
temperatures disagree with this single temperature fit at almost 
the 90% confidence level. However, the presence of the two 
point sources within the cooler sector and increased system- 
atic uncertainties in fitting the spectra of regions which are not 
azimuthally symmetric argues for caution in interpreting these 
results. We conclude that the uniformity of the gas temperature 
with azimuthal angle cannot be ruled out at the greater than 
90% confidence level. 

6. CENTRAL COOLING FLOW SPECTRA 

6.1. General Properties 

We also extracted the spectrum of the inner region of the 
cluster to study the properties of the cooling flow in more de- 
tail. The SIS and GIS spectra were extracted for the inner 3 
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arcmin (360 kpc) in radius of the cluster. This region just fits 
on the SOCl and S1C3 chips, which are the "best" chips on 
the two instruments and the ones normally used to determine 
the spectra of point sources. We note that our treatment of the 
central cooling flow spectra differed from that of the previous 
analysis of the global cluster spectrum (§ ^ or of the spatially 
resolved spectra (§ |5]) in several ways. First, we included the 
spectra from the two SIS instruments; these were excluded in 
the analysis of the global spectrum because the cluster extended 
beyond the field of view of the SIS. Second, we include the 
lower energy channels in the spectra, which were not included 
in the spatial analysis in § g. Because our interest here is in the 
soft X-ray emitting cooling flow, we extend the spectral analy- 
sis down to energies of 0.35 keV in the SIS and 0.55 keV in the 
GIS. 
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Fig. 7. — The temperature of the gas in A2029 at radii between 3 and 16 
arcmin as a function of the position angle around the center of the cluster. The 
position angle PA is measured from north to the east. The dashed horizontal 
line shows the best-fit single temperature model for the same spectra. The error 
bars are the 90% confidence regions. 



The energy-dependent PSF of ASCA can cause problems in 
determining the spectrum of any component of an extended X- 
ray source like a cluster. In § pi we corrected for these effects 
while simultaneously fitting simple models to the X-ray spec- 
trum in several regions of the the cluster Unfortunately, we 
will fit more complicated models to the spectrum of the cool- 
ing flow, and such a simultaneous fit is not practical with these 
models. Also, to study the cooling flow, it is important to ex- 
tend the spectral analysis down to lower energies where the PSF 
is not well-understood. We expect that the problems associated 
with the energy-dependent PSF will be much smaller for the 
cooling flow region than for the outer cluster regions, because 
of the large surface brightness peak associated with the cool- 
ing flow. To test this, we have determined the contribution of 
outer X-ray emission scattered into the 3 arcmin circle from 
outer regions, assuming the ROSAT PSPC image for the clus- 
ter and the best-fit spectral parameters as a function of position 
from § H In the GIS, we find that scattered light from outer 
regions contributes roughly 8% to the flux in the central region. 
The values are somewhat smaller for the SIS. There is also an 
energy-dependent loss of flux from the central region due to 
the PSF. During the spectral analysis, the program ASCAARF 
corrects for this effect. 

We also extracted the ROSAT PSPC spectrum from the same 



region. All of the spectra were binned to insure that there were 
at least 20 counts per spectral channel, so that the x^ test should 
give reliable results. 

Single temperature model fits to the GIS, SIS, and PSPC 
spectra are shown in rows 15-17 of Table |l[ Individually, these 
are probably acceptable fits. However, we could not find an 
acceptable common fit to all of the spectra. There are two in- 
consistencies between the fits to the different instruments. First, 
the absorbing column required by the SIS spectrum was about 
twice that found with the PSPC or the GIS. Second, the tem- 
peratures found by each of the instruments disagreed at a level 
beyond the uncertainties. It is interesting that the temperature 
increased with the hardness of the spectral response of the in- 
strument. This suggests that multiple temperature components 
are present, and that each instrument shows the gas which emits 
more strongly in its spectral band. 

The three instruments did give consistent values for the heavy 
element abundance, and the values for the central spectra are in 
good agreement with those determined for the entire cluster. 
This agrees with our previous conclusion in § 5.2 that the abun- 
dances are nearly constant within A2029. 

Because of the inconsistent temperatures found with the dif- 
ferent instruments and the other evidence for a cooling flow at 
the center of A2029 (§§ |i[ |, & |]l]), we also fit the single tem- 
perature plus cooling flow model to the spectra from each of 
the instruments. The best-fit spectral parameters are listed in 
rows 18-20 of Table |l]. The data, best-fit spectral models, and 
residuals for the three instruments are shown in Figures |[ ^ 
and |l^. In all three cases, the addition of a cooling flow sig- 
nificantly improved the fit. The abundances are all consistent 
with one another and with the single temperature models. The 
temperatures from the three instruments are now all marginally 
consistent within the uncertainties. The cooling rates are also 
all consistent, but the uncertainties are rather large. A cooling 
rate of about 250 Mq yr"' within the central 3 arcmin would be 
most consistent with the three individual spectral fits. 

Unfortunately, the required absorbing columns disagree even 
more strongly for the cooling flow model than they did in the 
single temperature model. When we tried to simultaneously fit 
the spectra from all three instruments with a single model, we 
were unable to find any acceptable fit. We allowed the nor- 
malization of the model to vary from instrument to instrument 
to account for any error in the overall calibration of the instru- 
ments, and only required that the model have the same spectral 
shape. Still, no acceptable common fit was found. The last row 
in Table [I] shows the best-fit common single temperature plus 
cooling flow model fit to all of the spectra. It is a rather poor fit. 

Examination of the residuals to such fits showed that the 
problem was primarily due to a difference in the shape of the 
spectrum and the fluxes at low photon energies (<1 keV) be- 
tween the ASCA SIS and the ROSAT PSPC. We tried a num- 
ber of different models for the origin of the soft X-ray absorp- 
tion (e.g., partial covering models, intrinsic absorption models, 
and models with intrinsic absorption applied only to the cool- 
ing flow). We also considered various alternative models for 
the cooling flow spectra (e.g., non-isobaric cooling or a finite 
limiting temperature for the cooling). None of these provided 
an acceptable common spectral fit or a significant improvement 
on the unacceptable common cooling flow fit. These problems 
suggest that there are significant calibration errors in the SIS 
data at low energies at present. Similar calibration problems 
and overestimated absorbing columns with the SIS have been 
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Abell 2029 
Cluster Center (r < 3 arcmin) CIS — Cooling Flow Mode 




energy (keV) 

Fig. 8. — The ASCA GIS2 and GIS3 X-ray spectra for the central cooling flow region for A2029 (r < 3') as a function of the 
measured photon energy are shown in the upper panel. The crosses give the data points with l-cr error bars, while the histogram 
is the best single temperature plus cooling flow model for the GIS spectra (row 18 in Table nl). The width of the data points or 
histogram steps is the width of the energy channels used to accumulate the data. The upper spectrum is for the GIS3, while the lower 
spectrum is for the GIS2. The GIS3 has larger energy bins because the bit error in memory during the observation reduced the energy 
resolution (§ 2. 1 ). The lower panel gives the residuals to the fit (in counts/sec/keV). 



found by other observers (e.g., Mukai 1995). These low energy 
SIS calibration uncertainties should not affect the previous re- 
sults (§§ H & pb, because low energy SIS data was not used 
there. 

6.2. Excess Soft X-ray Absorption? 

White et al. (1991) reanalyzed the Einstein Observatory 
Solid State Spectrometer (SSS) spectrum of the central 3 ar- 
cmin radius of A2029 (essentially the same region as used to 
extract the ASCA and ROSATPSPC spectra). White et al. found 
an excess absorption column of ANh = 15 x lO^*' cm"^ (if the 
excess absorber is modeled as Galactic material). In their fits 
to the spectrum, the absorber was treated as lying in the fore- 
ground of all of the X-ray emission. 

Such a large foreground excess column is completely incon- 
sistent with the ROSAT PSPC spectrum. The 90% upper limit 
on the foreground excess absorbing column in the inner 3 ar- 
cmin of A2029 from the PSPC spectrum is ANh < 2.05 x 10''* 
cm"^. Figure iri| shows that best-fit single temperature plus 
cooling flow model for the PSPC spectrum of the inner 3 ar- 
cmin radius of A2029, with the absorption fixed at the value 
from White et al. (1991). This is a terrible fit to the spectrum, 
with a x^ of 2417 for 23 d.o.f. The strong soft X-ray emission 
in the spectrum from the center of the A2029 cluster is incon- 
sistent with any large amount of excess foreground absorption. 
The White et al. value for the excess foreground absorption is 
also inconsistent with the ASCA SIS and GIS spectra of the cen- 
tral 3 arcmin of A2029 (rows 18 and 19 in Table Rb, albeit at a 



lower level of significance. 

If the excess absorption is strongly concentrated to the cen- 
ter of the cooling flow, then its effects would be diluted in the 
3 arcmin aperture used to collect the cooling flow spectrum. 
In order to provide a stronger constraint on excess absorption, 
we also extracted the ROSAT PSPC spectrum from the inner 1 
arcmin radius of the cooling flow. To provide a conservative up- 
per limit on any foreground excess absorption, we compared the 
absorption derived from this inner region with both the Galac- 
tic column of A^H = 3.07 x 10^° cm'^ from Stark et al. (1992), 
and the galactic absorption derived from the PSPC spectrum of 
the outer parts of the cluster (the annulus from 3-16 arcmin in 
radius). We adopted the larger value of the upper limit on the 
absorption for these two determinations of the Galactic column. 
The excess absorber was treated as a foreground screen, but at 
the redshift of the cluster (this increases slightly the limits on 
the excess column). This gave a 90% upper limit on the column 
of foreground excess absorption of 



ANniforeground) < 7.3 x 10 



19 -2 

cm . 



(3) 



The detection of excess absorption toward the central regions 
of cooling flow clusters of galaxies would suggest that the ab- 
sorber is physically located with the cooling flow region (White 
et al. 1991; Aflen et al. 1993; Irwin & Sarazin 1995). In this 
case, it is much more difficult to detect or constrain the excess 
absorption. To illustrate this, we consider a model in which the 
excess absorption is applied only to the cooling flow compo- 
nent of the spectrum, but not to the ambient cluster emission. 
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Abell 2029 
Cluster Center (r < 3 arcmin) SIS — Cooling Flow Mode 




energy (keV) 

Fig. 9. — The ASCA SISO and SISl X-ray spectra for the central cooUng flow region for A2029 (r < 3')- The notation is the same 
as for Figure pi The histogram is the best single temperature plus cooling flow model for the SIS spectra (row 19 in Table 11]). The 
upper spectrum is for the SISO, while the lower spectrum is for SlSl. 



If this model is fit to the PSPC spectra of the inner 3 arcmin of 
A2029, the best-fit value of the excess absorption is still zero. 
However, the 90% upper limit on the excess absorption is very 
large, 

ANHicooling flow) < 2.6 x 10^' cm"^ . (4) 

The reason for this very weak limit is that the value of the ex- 
cess absorption and the cooling rate are very strongly correlated 
in the fits. If one increases the amount of soft X-ray emission 
by increasing the cooling rate Mcoh one can increase the excess 
soft X-ray absorption /SNh if it is applied only to this cooling 
flow emission. 

Our conclusion is that the PSPC spectra (and, to a lesser ex- 
tent, the ASCA spectra) rule out the large foreground absorption 
claimed by White et al. (1991). However, an absorber confined 
to the cooling flow region is much more difficult to detect or 
constrain with PSPC spectra. In principle, ASCA SIS spectra 
could give a stronger constraint, but the inconsistency between 
the ASCA SIS and ROSAT PSPC spectra in A2029 suggests 
that there are still significant calibration uncertainties with the 
ASCA SIS spectra at low energies. 

For the purposes of determining the X-ray gas pressure in the 
inner region where the radio source is located (Huang & Sarazin 
1998), we also determined the best-fit PSPC temperature for the 
inner arcmin. This gave kT ■ 



:3.16!!Jj^i.eV. 



7. MASS PROFILE 



We used the information on the temperature and temperature 
distribution from the ASCA spectra to determine the gas and to- 
tal mass distribution in A2029. The first step in this process was 



to determine the X-ray surface brightness distribution, using the 
ROSAT PSPC image (Fig. nl). The X-ray count rate was deter- 
mined in concentric circular annuli, excluding regions around 
each of the point sources in the field. The surface brightness 
was corrected for background using the sa me r egions used for 
the background correction of the spectra (§ 2.2). The resulting 
X-ray surface brightness profile for the spectral band 0.5-2.0 
keV is shown in Figure n3. A sufficient number of counts were 
collected in each annulus to insure that gaussian statistics ap- 
plied. Possible faint emission is seen extending to a radius of 
at least 22 arcmin (2.64-h'^l^ Mpc), but the uncertainties in the 
surface brightness are very large, due to the low surface bright- 
ness and the effect of the instrumental support structure in this 
region. In the subsequent analysis (and in Fig. n2b, we consider 
only the surface brightness at radii r < 1 8 arcmin (2 . 1 6/15,' Mpc) 
where the uncertainties are more reasonable. 

We fitted the X-ray surface brightness in Figure O with the 
isothermal "beta" model 



h{r) = Io 



1 + 



-1 -3/3+1/2 



(5) 



The model was convolved with the PSPC PSF and compared to 
the data. When all of the data points in Figure n3 were included, 
no acceptable fit could be found. The best-fit model, shown as a 
dashed curve in Figure O, had x^ = 1074.5 for 27 d.o.f. As has 
generally been found with cooling flows clusters (e.g., Jones & 
Forman 1984), the beta model does not fit the sharply peaked 
central surface brightness within the cooling radius. Thus, we 
fit the surface brightness with the beta model, progressively 
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Abell 2029 
Cluster Center (r < 3 arcmin) PSPC — Cooling Flow Mode 
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Fig. 10. — The ROSAT PSPC X-ray spectrum for the central cooHng flow region for A2029 (r < 3')- The notation is the same as for 
Figure pi. The histogram is the best single temperature plus cooling flow model for the PSPC spectrum (row 20 in Table ll|). 



removing the interior points until an acceptable fit (xVd.o.f 
< 1.4) was found. This required that the data interior to 164" 
(0.33/zJ() Mpc) be eliminated. When these points were removed, 

108!?? arcsec and f3 = 0.64 ± 0.02 



an acceptable fit with rcore - ^^u-ie 
was found. Our value of (3 is in very good agreement with that 
given by David et al. (1995) of (3 = 0.65!U:Jj^. Our fit, which 
gave x^ = 18.9 for 16 d.o.f. (xVd.o.f = 1.182), is shown as the 
solid curve in Figure |l2[ Because the derived core radius is 
within the region of data excluded from the fit and where the 
beta-model doesn't fit, the value of the core radius and central 
surface brightness I,, are probably not given reliably by this fit. 
However, this model does appear to provide a reasonable fit for 
the surface brightness beyond about 160 arcsec. Because we 
are most interested in the mass at large radii, we will use this fit 
to describe the gas distribution. 

We adopt the linear fit to the radial temperature in the gas de- 
rived from the ASCA spectra (eq. H). Assuming the abundances 
of the heavy elements are 0.40 of the solar value, which was the 
best-fit abundance for the global cluster spectrum, we inverted 
the best-fit beta-model to determine the density of the gas pga., 
and its pressure P. The density of the gas was integrated over 
the interior volume to give the interior gas mass as a function of 
the radius, Mgasif)- The total gravitational mass was determined 
from the assumption that the gas is in hydrostatic equilibrium, 
which implies that 



Mror(r) = - 



dP 



Gpgas dr 



(6) 



The uncertainties in Mgas{r) and M,ot(r) were determined from 
the uncertainties and covariances in the coefficients of the fit to 



the temperature as a function of radius (eq. Hand Fig. H) and in 
the PSPC surface brightness fit (eq. || and Fig. Q. 

The resulting mass profiles are shown in Figure n3. At a 
radius of 16' (1 .92/150 Mpc, the largest region used for the tem- 
perature determination), the masses are Mgas = (2.52 ±0.77) x 



10^%^^ Mq and M,„, = (9.42 ±4.22) x lO'^^/jg^M©, and the 
ratio of masses is Mga/M,„, = (0.26 ± 0.14)/zgj^^. The fraction 
of the total cluster mass which is due to the intracluster gas is 
rising steeply with radius at the largest radii. A very similar 
result was found for A2029 by David et al. (1995) assuming 
an isothermal gas profile, and for Abell 2256 by Markevitch & 
Vikhhnin (1997) using an ASCA temperature profile. 



8. CONCLUSIONS 

We have analyzed the ASCA and ROSAT PSPC X-ray ob- 
servations of the rich cluster Abell 2029. The ASCA GIS and 
ROSAT PSPC are in good agreement on the global X-ray spec- 
trum, and give an average ambient gas temperature of 9.35!q45 
keV (including the effects of a cooling flow). If the gas temper- 
ature is assumed to be constant in the analysis of the spatially 
resolved ASCA SIS and GIS, the average temperature is found 
to be 8.6 keV. The iron abundance in the gas is 0.40 ±0.04 of 
the solar value. There is no significant evidence for any varia- 
tion in the abundance with position in the cluster 

The global X-ray spectra, central X-ray spectra, and ROSAT 
surface brightness all require a cooling flow at the cluster cen- 
ter The global X-ray spectrum implies that the total cooling 
rateis363!^^/!5^M0yr-i. 

The global X-ray spectra are consistent with the Galactic 
value for the soft X-ray absorption toward the cluster. The 
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Fig. 1 1. — The ROSAT PSPC X-ray spectrum for the central 3 arcmin radius of A2029 compared to the best-fit single temperature 
plus cooling flow model with the absorbing column set at the value from White et al. (1991). The notation is the same as in Figure nfl 



ROSAT PSPC spectra of the central regions of the cluster of 
completely inconsistent with the large value of foreground ex- 
cess absorption found by White et al. (1991) based on the Ein- 
stein SSS spectrum. The upper limit on excess foreground ab- 
sorption is 7 . 3 X 10''^ cm~^ . However, the spectra do not rule out 
a significant amount of intrinsic absorbing gas located within 
the cooling flow region. One concern with modeling the cool- 
ing flow spectrum and absorption in A2029 is that the ASCA 
SIS and ROSAT PSPC spectra do not appear to be consistent. 
This may indicate that there are calibration problems with the 
ASCA SIS spectra at low energies. 

The ASCA spectra of the cluster indicate that the gas temper- 
ature declines with radius. Similar declines are seen in most 
clusters studied with ASCA (Markevitch et al. 1997). 

The PSPC image shows that the cluster is elliptical, but is 
very regular and smooth. This agrees with previous analyses 
of the optical and X-ray distribution of the cluster (Slezak et 
al. 1994; David et al. 1995; Buote & Tsai 1996). We also find 
that there is no significant evidence for any irregularities in the 
temperature distribution in the cluster, as would be produced by 
a subcluster merger Structure in the X-ray surface brightness 
and particularly the gas temperature distribution has been inter- 



preted as evidence for mergers and other hydrodynamic activity 
in many other clusters (e.g., Henry & Briel 1995; Markevitch 
et al. 1996,1997). The lack of structure in the X-ray properties 
of A2029 suggests that the cluster is relaxed and that the gas is 
in hydrostatic equilibrium. 

We use the assumption of hydrostatic equilibrium to de- 
termine the gravitational mass of the cluster as a function of 
radius. Within a spherical radius of 16' (1.92/25^ Mpc), the 
total gravitational mass is M,o, = (9.42 ±4.22) x lO'^^/ijoM©, 
while the mass of gas is Mj,„ = (2.52 ±0.77) x lO^'^h-^MQ. 
The gas fraction of the cluster at this radius is Mga/M,o, = 
(0.26 ±0.14)/!5^, and the gas fraction is is increasing with ra- 
dius at the largest radii. Thus, A2029 is a particularly strong 
example of the so-called "baryon catastrophe" in clusters (e.g., 
David et al. 1995). 
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Fig. 12. — The ROSAT PSPC X-ray surface brightness profile for the spectral band 0.5-2.0 keV. The horizontal error bars show 
the width of the annulus used to collect the counts. The vertical error bars give the 90% confidence region. The surface brightness 
is plotted at the emission-weighted median radius for that annulus. The dashed curve is the unacceptable beta-model fit to the sur- 
face brightness including all of the data points. The solid curve is the acceptable beta-model fit, excluding the cooling flow regions 
(r < 164"). 
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Fig. 13. — The total gravitating mass and gas mass as a function of radius determined from the ASCA spectra and ROSAT PSPC 
surface brightness profile of A2029. The upper and lower solid curves give the total mass and gas mass, respectively. The dashed 
curves give the 90% confidence interval for each of these quantities. The dash-dot curve near the upper solid curve gives the total 
mass if the gas is assumed to be isothermal at the best-fit single temperature, rather than having a linear temperature variation with 
radius. The radial scale is given in arcsec at the bottom of the plot, and in kpc at the top. All values assume Hg = SQh^o km s"' Mpc~' . 
The masses may be inaccurate at the smallest radii because of the effect of the cooling flow. 



